HCN Survey of Normal Spiral, IR-luminous and Ultraluminous 

Galaxies 

Yu Gao 1,2 ' 3 and Philip M. Solomon 4 
ABSTRACT 

We report systematic HCN J=l-0 (and CO) observations of a sample of 
53 infrared (IR) and/or CO-bright and/or luminous galaxies, including seven 
ultraluminous infrared galaxies, nearly 20 luminous infrared galaxies, and more 
than a dozen of the nearest normal spiral galaxies. This is the largest and most 
sensitive HCN survey of galaxies to date. All galaxies observed so far follow 
the tight correlation between the IR luminosity Lir and the HCN luminosity 
L H cn initially proposed by Solomon, Downes, & Radford, which is detailed in a 
companion paper. We also address here the issue of HCN excitation. There is no 
particularly strong correlation between -Lhcn and the 12/zm luminosity; in fact, 
of all the four IRAS bands, the 12/xm luminosity has the weakest correlation 
with the HCN luminosity. There is also no evidence of stronger HCN emission 
or a higher ratio of HCN and CO luminosities Lhcn/^co f° r galaxies with 
excess 12/xm emission. This result implies that mid-IR radiative pumping, or 
populating, of the J=l level of HCN by a mid-IR vibrational transition is not 
important compared with the collisional excitation by dense molecular hydrogen. 
Furthermore, large velocity gradient calculations justify the use of HCN J=l-0 
emission as a tracer of high-density molecular gas (£3 x 10 4 /r cm -3 ) and give an 
estimate of the mass of dense molecular gas from HCN observations. Therefore, 
L H cn m ay be used as a measure of the total mass of dense molecular gas, and 
the luminosity ratio Lhcn/^co m ay indicate the fraction of molecular gas that 
is dense. 
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1. INTRODUCTION 

Though CO traces most of the molecular gas mass in galaxies, it does not necessarily 
trace active star-forming regions where the gas density is more than 10 times higher than 
the average. These regions are better traced by high dipole moment molecules like HCN, 
CS and HCO + . Consequently, these molecules have often been used as probes of physical 
conditions in giant molecular cloud (GMC) cores — the regions of active star formation in 
the Milky Way. Yet there have been few systematic studies of dense gas in GMCs in the 
disk of the Milky Way. Lee, Snell, & Dickman (1990) have found that for Milky Way disk 
GMCs, the ratio of CO/CS intensity is 156, much larger than for Galactic center clouds and 
an order of magnitude larger than for the archetypal ULIG Arp 220 (Solomon, Radford, 
& Downes 1990). Some dense clumps in some nearby GMCs have also been surveyed in 
CS (Lada, Bally, & Stark 1991; Plume et al. 1997) and in HCN (Heifer & Blitz 1997b). 
Although there are some recent attempts in Galactic plane HCN, CS, and CO observations 
(Heifer & Blitz 1997b; McQuinn et al. 2002 and references therein), only the central region 
of the Galaxy has been extensively mapped with dense molecular gas tracers (e.g., CS, Bally 
et al. 1987, 1988; HCN, Jackson et al. 1996; Lee 1996). Systematic and unbiased large-scale 
observations of the dense molecular gas in Galactic disk GMCs need continued study since 
the large-scale (kpc) distribution of dense gas in the Milky Way disk (cf. McQuinn et al. 
2002) is a basis for comparison with dense molecular gas observations in external galaxies 
(e.g., Paglione et al. 1998). 

HCN and CS are probably the most frequently observed interstellar molecules after 
CO. Because of their higher dipole moments (/i ~ 2.0 - 3.0 D), they require about 2 
orders of magnitude higher densities for collisional excitation than CO (/i ~ 0.1 D). HCN 
is one of the most abundant high dipole moment molecules and traces molecular gas at 
densities n(H 2 )~3 x 10 4 /r cm -3 (compared to densities at ^500 cm~ 3 traced by CO). 
HCN J=l-0 emission has been so far detected toward the nuclei of ~ 10 nearby galaxies 
(Nguyen-Q-Rieu, Nakai, & Jackson 1989; Nguyen-Q-Rieu et al. 1992; Henkel et al. 1990), 
10 relatively distant luminous and ultraluminous infrared galaxies (LIGs and ULIGs with 
IR luminosities 5 Ltr > 10 11 L and > 10 12 L , respectively) including three less luminous, 



5 The total IR luminosity (from 8 /im to 1000 ^m) is calculated according to the prescription in Sanders 
& Mirabel (1996) using measurements from all 4 IRAS flux bands. 
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"normal" spiral galaxies (Solomon, Downes, & Radford 1992), 12 nearby starburst and/or 
normal galaxies (Heifer & Blitz 1993), 10 interacting galaxies and infrared galaxy mergers 
(Aalto et al. 1995), 13 Seyfert galaxies (Curran, Aalto, & Booth 2000), and a few nearby 
CO-bright galaxies (Israel 1992; Sorai et al. 2002; Kuno et al. 2002). There are significant 
overlaps in the galaxy samples among these various observations, and the total number of 
galaxies with HCN detections is still ;S30, particularly the total HCN emission measured 
globally from the galaxies. A recent summary of HCN observations in centers of nearby 
galaxies can be found in Sorai et al. (2002) and Shibatsuka et al. (2003). 

In a few of the nearest galaxies, interferometer HCN maps have also been obtained in 
the nuclear regions (Carlstrom et al. 1988; 1990; Radford et al. 1991; Downes et al. 1992; 
Brouillet & Schilke 1993; Jackson et al. 1993; Heifer & Blitz 1995, 1997a; Tacconi et al. 
1994; Shen & Lo 1995; Paglione, Tosaki, & Jackson 1995b; Kohno et al. 1996; Kohno, 
Kawabe, & Vila-Vilaro 1999). CS emission, another dense molecular gas tracer, has been 
observed in the central regions of ~ 10 galaxies as well including interferometry imaging 
(Mauersberger & Henkel 1989; Mauersberger et al. 1989; Sage, Shore, & Solomon et al. 
1990; Heifer & Blitz 1993; Xie, Young, & Schloerb 1994; Paglione et al. 1995a; Peng et al. 
1996; Reynaud & Downes 1997; Wild & Eckart 2000). The presence of large amounts of 
dense molecular gas in the central ~ 1 kpc of galaxies has been well established (for an 
earlier review, see Mauersberger & Henkel 1993). 

The main goal of this paper is to present a systematic HCN survey of a large sample 
of normal spiral galaxies and relatively distant LIGs and ULIGs. This includes an attempt 
to measure the total HCN emission in nearby large spiral galaxies by mapping the HCN 
distribution along the major axes (Gao 1996, 1997). We examine correlations between 
HCN emission and various IRAS bands and discuss the excitation mechanisms of HCN. 
We here argue that HCN is primarily excited by collision with molecular hydrogen at high 
density and thus can be used as a tracer of dense molecular gas. We concentrate here on 
presenting observations of HCN J=l-0 only since it is almost always the strongest line 
among these high dipole moment molecules, and thus the easiest and most often observed 
dense gas tracer in external galaxies (Gao 1996). We have not only doubled the total 
number of galaxies with HCN J=l-0 detections in the literature for the last decade, but 
we have also extensively mapped the inner disk HCN emission in a few nearby galaxies in 
order to estimate the total HCN emission. We will explore in detail the extent and radial 
distribution of dense molecular gas in nearby galaxies in a future paper and present the 
detailed HCN mapping observations there (Y. Gao & P.M. Solomon in preparation). In a 
companion paper (Gao & Solomon 2004, Paper II), we fully discuss the interpretation and 
implications of our HCN survey. 
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We describe the sample selection and characteristics of our HCN survey in the next 
section. In § 3, we discuss our various observations and the observational strategies 
conducted at different telescopes. The survey results are presented in § 4, where we show all 
the HCN J=l-0 and CO spectra, tabulate the observational data, and summarize the global 
quantities for all galaxies. We have confirmed the tight correlation between the IR and 
HCN emission, initially proposed by Solomon et al. (1992) based on the HCN observations 
of a sample of 10 representative galaxies, in a large and statistically significant sample of 
more than 60 galaxies when both samples are combined. We have also checked for any 
significant dependence of the HCN emission upon the excess mid-IR emission as indicated 
by an IRAS flux ratio of the 12 and 100 /iin emission or the ratio of the 12 /im and total 
IR emission to address the possibly enhanced excitation of HCN J=l-0 emission by mid-IR 
radiative pumping through a vibrational transition. Moreover, the use of HCN as a tracer 
of high-density molecular gas is discussed using large velocity gradient (LVG) calculations. 
Finally, we summarize the main points of this HCN survey in § 5. 

2. THE HCN SURVEY SAMPLE 

Our goal is to survey a large sample of galaxies (> 50) with a wide range of IR 
luminosities where the total HCN J=l-0 emission from the whole inner disk or entire 
galaxy can be measured. Essentially, all galaxies with strong CO and IR emission were 
chosen for our HCN observations. Our HCN survey sample is mainly drawn from the CO 
observations of Solomon et al. (1997), Sanders, Scoville, & Soifer (1991), Tinney et al. 
(1990), Solomon & Sage (1988), and Young et al. (1989, 1995). Almost all galaxies with an 
antenna temperature measured by various telescopes to be greater than 100 mK (or 50 mK 
for LIGs) were included as candidate galaxies. These cutoffs were specifically used in order 
to select all truly CO-bright galaxies in the hope of detecting their HCN emission (which 
is typically 10-20 times weaker than CO) with a rms sensitivity at the level of ~ 1 mK. 
Most IR-bright galaxies with feo^ m > 50Jy, or /loo^m > 100Jy are included as well in our 
sample (Rice et al. 1988; Strauss et al. 1992). We thus selected ~ 10 nearest CO and/or 
IR bright large galaxies with cz^l,000 kms -1 and D 2 5~10', and a few tens of relatively 
distant yet still nearby "normal" spiral galaxies and LIGs. Moreover, most northern ULIGs 
in the local universe (cz < 20,000 kms" 1 ) were also added to enlarge our sample at this 
extremely high IR luminosity. 

We detected HCN J=l-0 emission from 52 galaxies, and only one galaxy was not 
detected. Table 1 lists all 53 sources, along with some basic properties of the galaxies 
useful for our observations and data analysis. The nearest large galaxies that require 
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mapping to estimate the total HCN emission are highlighted in boldface in Table 1, and 
the detailed observational results of these galaxies will be discussed separately (Y. Gao 
& P.M. Solomon 2004, in preparation). Seven ULIGs are marked with a star symbol in 
Table 1. Our HCN survey sample also contains essentially a statistically complete sample 
for northern (5 > —35°) galaxies with /ioo^m > 100 Jy, which includes about 30 of the 
brightest IRAS galaxies. 

For the purpose of the detailed statistical analysis to be presented in Paper II, we 
included almost all HCN data available in the literature where the total HCN emission 
from galaxies can be estimated to augment our sample. There are about a dozen mostly 
IR-luminous galaxies with total HCN emission measured, almost entirely from Solomon 
et al. (1992), but M51 and NGC 4945 are from Nguyen-Q-Rieu et al. (1992) and Henkel 
et al. (1990), respectively. Even with these supplements and those intentionally added 
ULIGs, the survey sample is not severely biased toward the high IR luminosity end. As 
shown in Figure la, although most of our sample galaxies are IR bright or the brightest in 
the local universe, more than half have IR luminosities less than 1O 11 L , and among them 
most have Ltr < 3 x 1O 1O L . There seems to be a slight bias, however, toward high CO 
luminosity galaxies (Fig. lb), i.e., more CO-luminous or molecular gas-rich galaxies have 
been selected. This bias mainly arises because most normal spiral galaxies selected are 
CO bright and quite rich in molecular gas, and essentially all LIGs and ULIGs have high 
CO luminosities. Nevertheless, the galaxy number distribution in HCN luminosity appears 
to be rather unbiased with comparable numbers of galaxies with both low and high HCN 
luminosities, and the entire HCN luminosity distribution covers a large range of almost 3 
orders of magnitude (Fig. lc). 

3. OBSERVATIONS 

We used the former NRAO 6 12m telescope at Kitt Peak and the IRAM 30m telescope 
at Pico Veleta near Granada, Spain, for most of our HCN observations. The FCRAO 7 14m 
telescope with the 15 element focal-plane array receiver, QUARRY (a configuration of three 
rows of five beams with a row spacing of 100" and a beam spacing of 50"), was initially used 



6 The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under 
cooperative agreement with the National Science Foundation. 

7 The Five College Radio Astronomy Observatory is operated with support from the National Science 
Foundation, and with the permission of the Metropolitan District Commission, Commonwealth of 
Massachusetts. 
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for mapping some of the nearest strong sources. We conducted several observing runs at 
each telescope from 1993 April to 1995 December. Additional observations were conducted 
remotely at the NRAO 12m telescope in 1997 May and June to confirm some weak 
detections and to further improve the signal-to-noise ratios for some HCN spectra. The 
system temperatures at the redshifted HCN line frequencies were typically 220, 250, and 500 
K on average on the antenna temperature scales of T^, T\, and T\ for the NRAO, IRAM 
and FCRAO telescopes, respectively. We converted the observed antenna temperature to 
the main-beam temperature scale (T m b) for all our observations at different telescopes. 
Table 2 lists various telescope beam efficiencies measured at different telescopes at both 
the 2.6 and 3.4 mm wavebands. At ~ 3 mm, we adopted the continuum sensitivity (the 
temperature-to-flux conversion for a point source) on the T m b scale as S/T m ^ = 4.8, 29.0, 
and 23.5 Jy KT 1 for the IRAM 30m, NRAO 12m, and FCRAO 14m telescopes, respectively. 

At the IRAM 30m, we used two 512 x 1 MHz filter banks and the autocorrelator in its 
widest bandwidth configuration as back ends. These back ends were connected to three SIS 
receivers (3mm, 2mm, 1mm) that were usually used to simultaneously observe HCN J=l-0, 
CS (3-2), and CO J=2-l lines respectively. The telescope beam diameters (FWHM) for CO 
J=l-0 and HCN J=l-0 were 22" and 27", respectively. The NRAO 12m was equipped with 
dual-polarization 3mm SIS receivers. We connected the dual receivers to a pair of 256 x 2 
MHz filter banks as well as the two 600 MHz wide hybrid spectrometers. The telescope 
beam diameters (FWHM) for CO J=l-0 and HCN J=l-0 were 55" and 72", respectively. 
At the FCRAO 14m, the corresponding back ends for the 15 cooled Schottky diode mixers 
of the QUARRY array were 15 spectrometers of 64 x 5 MHz bandwidth. The telescope 
beam spacing of 50" within each row of QUARRY is about the beam size (FWHM) of the 
14m telescope at 3 mm. 

The IRAM 30m telescope used a wobbling secondary with a beam throw of ±4' to 
achieve flat baselines. Similarly, a nutating subreflector with a chop rate of ~ 1.25 Hz 
was used at the NRAO 12m together with a position switch (the so-called beam-switch 
plus position-switch, or "BSP mode"). The beam throws used were in the range of ±2' 
to ±3' (the larger beam throw was used for nearest large galaxies with D 2 5~10'). Thus, 
the separations between the OFF reference positions and the ON source position are in 
the range of ±4' to ±6'. For a couple of the largest galaxies (_D25~20') observed at the 
12m, we employed position switching in azimuth with offsets comparable to the expected 
CO source sizes to ensure that the reference positions are completely outside of the sources 
observed. At the FCRAO 14m, we used double-position switching in azimuth with offsets 
comparable to the diameters of galaxies. The QUARRY receivers of the FCRAO 14m 
telescope were efficient in mapping the CO emission in the inner disks of nearby galaxies, 
but the sensitivity achievable limited our extensive mapping of the weak HCN emission 
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outside the nuclear regions. No additional beams except for the central beam and sometimes 
the adjacent two inner beams were detected in HCN with the 14m. Therefore, most of the 
HCN mapping was eventually conducted with the 12m and 30m telescopes after the initial 
trial observing runs with the QUARRY at the FCRAO 14m. 

We mainly used the IRAM 30m to observe most (rather distant) ULIGs in the sample 
and to map some nearby starburst galaxies of mostly smaller optical diameters, given the 
smaller telescope beam compared with the source size and extent for achieving effective 
mapping. Essentially all other observations were conducted with the NRAO 12m so that 
one beam measurement covered all emission from relatively distant galaxies including some 
interacting/merging galaxy pairs. Almost all nearby spiral galaxies were mapped along 
the major axes with the NRAO 12m. These 12m observations included both the HCN 
observations initially attempted with the FCRAO 14m QUARRY array, and in several 
cases, the major axis mapping observations from the IRAM 30m in order to properly 
compare the total HCN emission estimated with different telescopes (see Table 3). This 
is essential to better characterize the radial HCN distribution in galaxies (Y. Gao & P.M. 
Solomon 2004, in preparation). 

All observed positions were chosen by checking and comparing the positions of the 
peak radio continuum emission as determined from the VLA snapshot survey (Condon 
et al. 1990) or the best positions (with uncertainties ^10") then listed in the NED. 8 
The observed spectra from all telescopes (except for the FCRAO 14m) were smoothed 
to ~8 MHz resolution, or ~ 20 kms -1 and ~ 27 kms~ x at 2.6 and 3.4mm, respectively. 
Occasionally, in several of the weakest HCN sources, the spectra were smoothed to 16 MHz 
to boost the signal-to-noise ratio. Linear baselines were removed from the spectra after 
emission line windows were chosen to estimate the integrated line intensities. We converted 
the line temperatures to the main-beam brightness temperature scale T mb (Tables 2 and 3). 
In many cases, CO lines were also observed from the same telescope for comparison and 
for determination of the HCN line emission windows. All HCN observations of galaxies so 
far have shown that the significantly detected HCN line profiles have the same line widths 
as the CO lines (e.g., Solomon et al. 1992; Heifer & Blitz 1993; Gao 1996). With the 
CO spectra obtained from the same positions or available in the literature from the same 
telescope, the integrated line intensities of the weakly detected HCN lines can be more 
accurately estimated. 

Calibrations at all three telescopes were done using the standard chopper wheel method 



8 Thc NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space Administration. 
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(Kutner & Ulich 1981). The absolute antenna temperature scale was checked by either 
observing the standard calibration sources or by the measured CO line strengths of the 
targets (mostly for the nearest strong sources) as compared with the standard spectra at 
the telescopes and the known CO spectra available in the literature. Pointing was checked 
by observing quasars or planets about every 2 — 3 hours, with typical errors of ~ 3" for 
the IRAM 30m, and ~ 5" — 10" for the NRAO 12m and the FCRAO 14m. Focus was 
also monitored a few times per day mostly using planets. More detailed discussions of the 
observations and the antenna temperature scales at different telescopes can be found in 
Gao (1996, 1997). 



4. RESULTS AND DISCUSSION 
4.1. Spectra 

We present all HCN spectra in Figs. 2 and 3 for the IRAM 30m and the NRAO 12m 
telescopes, respectively, including a few weak (~ 3a) detections. For most galaxies, we also 
show the CO spectra taken at the same position from the same telescope for comparison. 
A few CO spectra in Figure 2, however, are CO J=2-l lines, which have much smaller 
beam sizes compared with that of the HCN J=l-0 spectra. In some cases, the velocity 
range of CO emission, obtained from the published CO spectra in the literature with 
the same telescopes or comparable telescope beam sizes, is indicated instead. NGC 253, 
NGC 2903, M82, and NGC 6946 were mapped by both the IRAM 30m and the NRAO 12m 
and therefore appear in both Figs. 2 and 3. Only the central HCN spectra of those nearby 
galaxies highlighted in Table 1 are displayed here; the detailed HCN mapping will be shown 
in a future paper (Y. Gao & P.M. Solomon 2004, in preparation). 

We detected HCN emission in essentially all of the sources (except NGC 1055) listed in 
Table 1, though a few galaxies were only marginally detected. The observed line intensities 
(of both HCN and CO if applicable) are listed in Table 3 together with other measured 
quantities useful for analysis. The derived global properties of the line luminosities and 
their ratios using methods discussed below (§4.2) are listed in Table 4. 

Fully sampled mapping is essential in order to estimate the total HCN emission in large 
spiral galaxies (Gao 1996, 1997; Y. Gao & P.M. Solomon 2004, in preparation) because a 
substantial fraction of the total HCN emission is distributed outside the ~ 1 kpc nuclear 
regions (which are covered by the central beams in these observations). Therefore, we did 
not add more nearby galaxies with HCN detections from their centers from the literature 
into our HCN sample. Those data only measure lower limits to the HCN luminosity since 
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only the central beam detections were made for these nearby large spiral galaxies, and 
the HCN emission outside the central ~ 1 kpc nuclear regions is unknown. Notice that 
about 15% of the galaxies in our HCN survey sample also have either upper or lower 
limits in the total HCN emission or luminosity (Tables 3 & 4). The upper limits in the 
HCN measurements correspond to weak detections (~ 3c) or non-detections (£j2a, only 
for NGC 1055) or to nearby galaxies with good detections in the center but with only 
weak detections or non-detections at the off-center positions. The lower limits in the HCN 
luminosity are for a few nearby galaxies with good HCN detections from the galaxy centers, 
but no systematic mapping or detections at off-center positions (NGC 2903 is exceptional 
since it was mapped with the IRAM 30m, yet more sensitive and extensive mapping is still 
required to estimate the total HCN emission). This is similar to the case of more than 
a dozen nearby galaxies with only central beam HCN observations in the literature (thus 
lower limits to the total HCN emission). We include our own limits, but not the HCN limits 
in the literature in the final HCN sample. 



4.2. HCN and CO Line Luminosities 

For sources much smaller than the telescope beam, the HCN luminosity is calculated 
using 

L'hcn = j Tb(HCN) dV dQ s d A 2 » vr/(4/n2) 6 mh 2 J HC n d L 2 (1 + z)^ (1) 

(see Solomon et al. 1997 for details, 11 with a unit of K kms^ 1 pc 2 is introduced to 
distinguish it from normal L, which has proper power units), where Q s is the solid angle 
subtended by the source area, 

/hcn = / T mh (ECN)dV = j T b (HCN) dV dQ s /[(l + z)Q s * h ] (2) 

is the observed integrated line intensity (as converted to the beam-diluted main- 
beam temperature scale), T b (HCN) is the intrinsic source brightness temperature, 
fi s *b ~ 7r /(41n2)^ m b is the solid angle of the source convolved with the diameter 
(FWHM) of the telescope Gaussian beam 9 mb , and d^ and d A are the luminosity 
distance 9 and the angular size distance to the source, respectively, related by 
d L = d A (l + zf = 7994.5[1 + z-(l + zf/ 2 } Mpc. 

For large galaxies where HCN was observed at several positions along the major axis, 



9 i?o=75kms 1 Mpc 1 and qo — 0.5 were used in this paper 
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the total integrated HCN flux over the source, 

^hcn = / T mh {KCN)dVdn a . h , (3) 

is estimated following the empirical formula adopted in Solomon & Sage (1988; also 
Sage 1987) to calculate the total CO emission. The total HCN emission is obtained by 
summing all the measured line fluxes along the major axes with weights determined from 
considerations of # m b, the (source) disk size, D 2 $, and the disk inclination, and assumption 
of an axisymmetric source distribution. Thus, for nearby galaxies, 

^hcn = / T mb (HCN) dV dQ s * h d A 2 (l + z) = Fhcn d L 2 (1 + z)~ 3 . (4) 

The CO luminosity is calculated similarly with the same formula listed above (switching 
HCN with CO). For brevity, we hereafter use the symbol L rather than V for the line 
luminosity in unit of K kms -1 pc 2 for both CO and HCN emission. 

For those nearby galaxies that we observed with different telescopes, the agreement 
in the estimated total HCN fluxes is rather good (Table 3), given the weak HCN lines 
detected outside the nuclear regions in most cases. For all positions observed in HCN, the 
permitted velocity ranges are available either from our own, or published, CO observations. 
As discussed in § 3, the weak detections or upper limits for the HCN line intensities can be 
better estimated when the HCN line emission windows are known/assumed from the CO 
lines, and thus the HCN results are rather robust. The uncertainties in the total HCN (and 
occasionally CO) luminosities are up to ~ 50% including the calibration errors of different 
telescopes, but in most cases only ~ 30% or better. The largest uncertainties are only 
present in a few galaxies with the marginal detections (~ 3a) and in some nearby large 
galaxies with very limited points observed where extensive mapping along their major axes 
is still required for more accurate estimates of the total HCN emission. 

4.3. Correlation between HCN and IR Emission 

Figure 4 shows a tight linear correlation between Lrcn and Lir over 3 orders of 
magnitude in IR luminosity. The originally proposed correlation, based on 10 representative 
galaxies (Solomon et al. 1992), is also noted for comparison. We recalculated Lhcn, 
Leo (see $4.2) based on the measured line intensities Ihcn and Too listed in Table 1 of 
Solomon et al. (1992). Obviously, the new HCN data from a sample more than 5 times 
larger closely follow the trend established by Solomon et al. (1992) and have about the 
same scatters. A least-squares fit including published data (almost entirely from Solomon 
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et al. 1992), yet excluding our HCN limits, yields a slope of 1.00 ± 0.05, with a correlation 
coefficient R=0.94 (R 2 = 0.88). The correlation implies logL m = 1.00/ogL H cN + 2-9. 
When we include all galaxies in the fit and take the limits given in Table 4 as the HCN 
luminosities, the linear relation remains almost the same with logLm = 0.97/ogLHCN + 3.1 
and the same correlation coefficient. 

We have therefore confirmed the tight correlation between IR and HCN luminosities in 
a large, statistically significant sample. LIGs/ULIGs have the highest Lhcn and highest 
ratio of Lhcn/^cO) and many of them have Lhcn even larger than Leo of the Milky Way, 
which implies a relatively high fraction of dense molecular gas as well as an extremely 
large reservoir of dense molecular gas (discussed further in § 4.5). The detailed discussion 
of the tight correlation between L IR and Lhcn and other relationships, as well as their 
implications, will be fully explored in Paper II. In this paper, we would like to show that 
HCN emission is indeed a good tracer of dense molecular gas (in the next two sections). 
Therefore, the tight correlation between far-IR and HCN suggests that active star formation 
is strongly related to the mass of dense molecular gas and the main source of the strong IR 
emission in LIGs/ULIGs is mostly from active star formation. 

In Figure 5, we have divided the sample into two subsamples: LIGs and ULIGs with 
L\r > 1O 11 L are plotted as filled circles, and less luminous galaxies are shown as open 
circles. Apparently, all LIGs/ULIGs have Lhcn larger than 10 8 K kms -1 pc 2 , whereas 
the less luminous normal spiral galaxies have Lhcn~10 8 K kms^ 1 pc 2 . We also note 
that several ULIGs have Lhcn/Lcjo as high as ~ 0.25, whereas most normal galaxies 
have Lhcn/^co nearly 10 times smaller. Other systematic differences between these two 
subsamples can also be easily seen in various other correlations presented in Paper II. 

4.4. Excitation of HCN 

We here try to empirically demonstrate that the radiative pumping, or populating, 
of the J=l level of HCN through a vibrational transition at 14 /im wavelength, does 
not enhance significantly the excitation of the HCN J=l-0 emission. Since the mid-IR 
vibrational transition of HCN falls within the IRAS 12 fim bandpass and there are no 
other direct measurements of the 14 /zm emission available, we use the 12 /im flux as an 
approximation. We show in Figure 5 two correlation plots that involve the IRAS 12 /im 
emission to illustrate that the enhancement in the excitation of HCN through a vibrational 
transition, by this potential mid-IR radiative pumping, as postulated by Aalto et al. (1995), 
is not important. Previous HCN modeling also found that the radiative pumping from the 
dust continuum emission has little effect on the excitation of the rotational transitions in 
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the vibrational ground state of HCN (e.g., Stutzki et al. 1988; Paglione, Jackson, & Ishizuki 
1997). 

Although there is still a good correlation between the 12 /im and HCN emission (a 
squared correlation coefficient of R 2 = 0.75), the relationship is not particularly strong 
as compared with the much tighter correlation between IR and HCN emission (Fig. 4, 
R 2 = 0.88). In fact, the correlation between the 12/im and HCN emission appears to be 
the weakest among the correlations between the four IRAS waveband and HCN emission. 
Both the 60 and 100/im luminosities correlate much better with the HCN luminosity. For 
example, the correlation between 100/im and HCN emission (R 2 = 0.89) is significantly 
better than that between 12/im and HCN emission (R 2 = 0.75). A three-parameter fit of 
HCN emission as a function of both 12/im and 100/im emission clearly shows that most of 
the contribution to the fit is from the 100/im emission alone, whereas there is little from 
the 12/im emission. Should the mid-IR pumping be an important effect in the excitation of 
HCN in the vibrational ground state of HCN, a tighter correlation between the 12/im and 
HCN luminosities than any correlations between the HCN and 60/iin, 100/im, and total IR 
luminosities would be expected. We observe the contrary. 

We can also use the luminosity ratio, rather than the single luminosity alone, to 
normalize the correlations in order to eliminate the distance and galaxy size dependencies 
inherent in the luminosity correlations (more examples are in Paper II). There is only a 
weak correlation between Li 2fim _/ Lqq and Lhcn/-^co (R 2 — 0.29), whereas the correlation 
between £ioo/im/-kco and Lhcn/^co * s s ^ m fairly strong (R 2 = 0.56). Clearly, these results 
suggest that the mid-IR emission is only loosely correlated with the HCN emission, and 
the correlation between the far-IR emission and the HCN emission is significantly tighter 
in comparison. This essentially rules out the possibility that the mid-IR pumping in the 
vibrational ground state of HCN could be one of the significant contributors to the HCN 
J=l-0 emissivity. 

To strengthen these results, we compare the HCN luminosity and the ratio of HCN and 
CO luminosities (L H cn/-^co) with the excess 12 /im emission estimated by either the ratio 
of 12/im and total IR or the ratio of 12/im and 100/im (or 12/im and 60/iin) luminosities. 
We plot the luminosity ratio of L^m/Li^ versus Lhcn in Figure 5a and the luminosity 
ratios of Lhcn/^co versus £i 2Mm /Lioo Mm i n Figure 5b. Obviously, no correlations are 
observed in these plots and there is no sign of increased HCN emission or an increased 
Lhcn/-^co ratio for galaxies with a 12 /im emission excess. Similar plots for relationships 
between the luminosity ratios of Lhcn/^co and L^^m/^m or between Lhcn/^co and 
-^i2^m /-^60/im show no correlations either. Therefore, there is little evidence for a direct 
connection between the 12 /im and HCN emission, even for galaxies with excessive 12 /im 
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emission (Fig. 5). The contribution of mid-IR radiative pumping to the HCN excitation 
must be negligible when it is compared to the collisional excitation by dense molecular 
hydrogen (see § 4.5). 

Gao & Solomon (Paper II) show that Lhcn can be quite accurately predicted from 
the IR (or far-IR) and CO luminosities. In this three-parameter model fit, L IR appears to 
be slightly more important, but L co is also a significant predictor for L H cN[the HCN(L IR , 
Leo) model (see § 3.3.1 and Table 3 in Paper II)]. Here we can add L 12Atm as an additional 
parameter into the model fit for Lhcn to check whether Li2 M m is a non-negligible predictor 
for Lhcn as well. This could test the significance of the mid-IR emission contribution 
to HCN as compared with that of the far-IR and CO emission. We find that the model 
fit changes only marginally with no improvement in the correlation coefficient or the rms 
scatter of the fit. In fact, no difference is noticeable in the model fit whether is 
added or not. Therefore, even though there exists a correlation between 12/zm and HCN 
emission, the correlation is significantly worse than the correlations between the HCN and 
IR emission of three other IRAS bands, and the contribution to HCN from 12 /zm emission 
is negligible. We conclude that there is no clear evidence to claim that the mid-IR pumping 
can significantly enhance the HCN emission. 

For the archetypal AGN/starburst hybrid galaxy NGC 1068, which has the strongest 
mid-IR emission in our HCN sample, the Lhcn/Lco ra tio is also one of the highest. In 
particular, a very high HCN/CO ratio (^,0.3) within ~ 100 pc nuclear region is observed 
(Jackson et al. 1993; Tacconi et al. 1994; Heifer & Blitz 1995) in the interferometric maps, 
which appears to be the mid-IR compact core around AGN (Le Floc'h et al. 2001). Is this 
an example that the mid-IR pumping has significantly enhanced the global HCN emission? 
Apparently not, since the innermost ;$100 pc nuclear region around the AGN contributes 
little to the average HCN/CO ratios in the circumnuclear (~0.5-l kpc) starburst ring 
where most of molecular gas is located. Also in NGC 1068, the ratio /i2^m//iooMm = 0.16 
is much higher than in any other galaxy in the sample. Yet, it is still not certain whether 
the extreme excess of the mid-IR emission contributes significantly to the excitation of 
HCN. This is because galaxies with the lowest L 12Mm /Lm ratios can have even larger HCN 
luminosities than NGC 1068 (Fig. 5a) and galaxies with the lowest fufim/ /loo^m ratios 
appear to have similar or even higher Lhcn/Lco ratios than NGC 1068 (Fig. 5b). Thus, 
globally the high Lhcn / -^co r& tio and high Lhcn in galaxies may have nothing to do with 
the mid-IR emission excess. Again, there is no indication that the strong mid-IR emission 
is linked to the high ratio of HCN/CO or large HCN luminosity even for galaxies with the 
most extreme mid-IR emission excesses. 
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4.5. Mass of Dense Molecular Gas Traced by HCN 

The mass of molecular gas in galaxies is typically derived from the CO luminosity 
using the "standard" CO-to-H 2 conversion factor 

M(H 2 ) = 4.78 L co M (K kms' 1 pc 2 )" 1 (5) 

(e.g., Solomon & Barrett 1991; Young & Scoville 1991 and references therein), estimated 
from the Galactic disk GMCs in the Milky Way. Although this conversion corresponds 
to iV(Lf 2 )//co = 3 x 10 20 cm -2 (K kms -1 ) -1 , rather than the recent calibration of 
N(H 2 )/I co = 1.7 x 10 20 cm' 2 (K kms" 1 )- 1 (Dame et al. 2001; Hunter et al. 1997), 
we use it simply for easy comparison with previous CO studies of galaxies. As discussed 
in Downes, Solomon, & Radford (1993), Solomon et al. (1997), and Downes & Solomon 
(1998), however, we caution that it might overestimate the molecular gas mass in ULIGs 
by factors of 3-5. Unlike normal spiral galaxies, where a large fraction of CO emission 
comes mainly from GMCs distributed in the inner disk, the CO distribution in ULIGs is 
extremely concentrated in the nuclear regions and may fill the entire volume since even 
the intercloud medium may be molecular (Solomon et al. 1997; Downes & Solomon 1998). 
High-resolution CO imaging of ULIGs has indeed shown extremely high concentration 
of CO emission in the central few hundred parsecs (e.g., Downes & Solomon 1998). For 
ULIGs, the conversion factor is thus lower than for normal spiral galaxies. Nevertheless, we 
still adopt the standard CO-to-H 2 gas mass conversion (eq. [5]) for convenience since 90% 
of the galaxies in our sample are not ULIGs. The ratio of L IR /M(H 2 ), often referred to as 
the star formation efficiency, is also listed in Table 4. 

Observations of H 13 CN (Nguyen-Q-Rieu et al. 1992; Wild et al. 1992) indicate that 
extragalactic HCN emission is also optically thick, just like CO emission. The conversion 
factor between HCN and the mass of the dense molecular gas, however, is poorly constrained 
as there is no direct calibration from GMCs or GMC core regions in the disk of the Milky 
Way. From observations of higher transition CS emission, Plume et al. (1997) find a 
typical GMC core has a diameter of 1 pc and mass of 3800M Q . Then the rms velocity 
of a virialized cloud with these parameters is ~ 5 kms -1 , consistent with the observed 
FWHM line width. We have carried out large velocity gradient (LVG) calculations that 
include radiative trapping for HCN J=l-0 to determine the mass of dense molecular gas 
from the measured HCN brightness temperatures (e.g., Kwan & Scoville 1975), using the 
collision rates given by Green & Thaddeus (1974). For an HCN abundance (relative to H 2 ) 
of 2 x 10~ 8 (e.g., Irvine, Goldsmith, & Hjalmarson 1987; Lahuis & van Dishoeck 2000) and 
a velocity gradient of ~ 5km s -1 pc -1 , this gives [HCN] /(dv/dr) ~4x 10 _9 (kms _1 pc -1 ) -1 . 
For a kinetic temperature of 20-50 K and the intrinsic HCN line brightness temperature 
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T b > 5 K, this gives N(H 2 )/I UCN & 2 x 10 21 cm" 2 (K kms" 1 )- 1 or 

M dense (H 2 ) £ 25L HC n M (K kms" 1 pc 2 )" 1 (6) 

for H 2 at a density ^10 4 /r cm -3 , where M densc (H 2 ) means the dense molecular gas as 
traced by HCN emission, and r£l is the optical depth of the HCN J=l-0 line. 

For higher intrinsic HCN line brightness temperature T b > 15 K, higher gas density 
n(H 2 )^3 x 10 4 /r cm -3 , and the same abundance and velocity gradient, however, LVG 
results give N(H 2 )/I U cn £ 1.3 x 10 21 cm" 2 (K kms" 1 )" 1 , or 

M dense (H 2 ) £ 15L HC N M (K kms" 1 pc 2 )" 1 . (7) 

Alternatively, for a virialized cloud core of an average density <n(H 2 )>= 3 x 10 4 cm -3 
and Tb ~ 35 K (e.g., Radford, Solomon, & Downes 1991), 

/XT w 1 / 2 

M densc (H 2 ) « 2.1-^^ Lhcn ~ 10 Lhcn M (K kms" 1 pc 2 )" 1 , (8) 

almost same as that obtained from LVG approximation (eq. [7]). 

The LVG calculations above give perhaps an upper limit for the conversion factor 
between HCN and the mass of dense molecular gas (cf. Mauersberger & Henkel 1993) 
since the much lower intrinsic HCN line brightness temperature T b > 5 or 15 K was used 
besides other approximations in the LVG estimates. Higher Tb indeed tends to result in a 
smaller conversion factor allowed by LVG (eq. [7]). HCN hot cores in regions of massive 
star formation usually have extremely high temperature up to hundreds of degrees (e.g., 
Boonman et al. 2001). Thus, it is likely that a much smaller conversion factor holds 
under these extreme conditions. Moreover, there are some other mechanisms possible 
to help excite HCN molecule, e.g., excess HCN abundance (Bergin, Snell, & Goldsmith 
1996; Lahuis & van Dishoeck 2000) resulting from the active star-forming environment 
and high supernovae events, and collisions with electrons. Although these are most 
likely insignificant (e.g., mid-IR radiative pumping discussed in previous section) and 
the dominant mechanism should still be collisions with molecular hydrogen, they could 
work collectively to lower the conversion factor given in equations (6) and (7). Other 
constraints such as the total molecular gas mass estimated from the CO observations, the 
higher frequency CO line observations, and the dynamical mass estimations (e.g., Solomon 
et al. 1997; Downes & Solomon 1998), all seem to point to a smaller conversion factor, 
particularly for LIGs/ULIGs. Thus, a lower conversion factor given in equation (8) is more 
likely. 
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In fact, since the luminosity ratio Lhcn/^co i n ULIGs can be as large as ~ 25% (e.g., 
Mrk 231, Mrk 273 and IRAS 17208-0014), the true conversion factor between HCN and 
the mass of dense molecular gas should obviously be much smaller than the upper limit 
allowed by the LVG calculation (eq. [7]) in order to reconcile with the total molecular gas 
mass estimated from the standard CO-to-H2 conversion (eq. [5]). Particularly for ULIGs, 
even equation (8) might still overestimate the dense molecular gas mass since T b could be 
indeed very high (£50 K) and the standard CO-to-H 2 conversion of equation (5) should be 
reduced by a factor of 5 in ULIGs (Downes & Solomon 1998). The main difference between 
normal spiral galaxies and ULIGs in the HCN-to-H 2 conversion, however, is the drastically 
different molecular line brightness temperatures. While T b ~ 10 K is valid for most normal 
spiral galaxies, a factor of 5 higher T b ~ 50 K was deduced for ULIGs undergoing extreme 
starbursts (Downes & Solomon 1998), which argues for a much lower HCN-to-H 2 conversion 
factor than those given in the above equations. 

Nevertheless, this also indicates that ULIGs with the highest Lhcn/^co ratio are 
likely to have most (£50%) of their molecular gas at high density (£3 x 10 4 /r cm -3 ) 
and that Lhcn/^co ratio can indeed be used as an indicator of the fraction of dense 
molecular gas in galaxies. On the other hand, Lhcn/^co ratio also indicates the average 
molecular gas density in galaxies. Given the large uncertainty of the HCN-to-H2 conversion, 
the same as is true for CO-to-H 2 conversion, we suggest that a conversion factor of 
M densc (H 2 ) / L UCN ~ 10 M Q (K kms -1 pc 2 ) -1 (eq. [8]) is the more reasonable one to use, 
and we adopt this relation to estimate the total dense molecular gas mass in galaxies as 
tabulated for L m / M dense (H 2 ) in Table 4. Accurate determination of the conversion factor 
between HCN and the mass of dense molecular gas needs further extensive studies including 
higher frequency HCN observations (e.g., Jackson et al. 1995; Paglione et al. 1997) and 
more detailed modeling. 



5. CONCLUSIONS 

We present systematic HCN J=l-0 observations (complemented with CO observations) 
of a sample of 53 IR/CO-bright and/or luminous galaxies, which is the largest HCN sample 
and most sensitive HCN survey of external galaxies so far. Ultraluminous infrared galaxies 
have the highest HCN luminosities, usually several times larger than the CO luminosity of 
our own Galaxy. Many luminous infrared galaxies have HCN luminosities comparable to 
the CO luminosity of the Galaxy. The ratio of HCN and CO luminosities can be as large as 
~ 25% in ultraluminous infrared galaxies, nearly an order of magnitude higher than most 
normal spiral galaxies. 
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We compared the HCN line emission with the far-IR emission (an indicator of the rate 
of high-mass OB star formation). All galaxies surveyed follow the tight correlation between 
IR and HCN luminosities initially proposed by Solomon et al. (1992). Since stars form only 
in dense clouds and the HCN-emitting gas (i.e., the dense molecular gas) and IR emission 
appear closely and physically related, then active star formation is the dominant source of 
IR emission in luminous and ultraluminous galaxies. This is fully discussed in Paper II. 

There is no particularly strong correlation between the 12 /im and HCN luminosities. 
The correlations between 100 /im (and 60 //m) and HCN emission are much better correlated 
in comparison. Galaxies with excess 12 /im emission do not show stronger HCN emission or 
higher HCN/CO luminosity ratios. Thus, mid-IR radiative pumping of HCN excitation is 
of little significance when compared with the collisional excitation by dense molecular gas. 

We also discussed the use of HCN J=l-0 emission as a tracer of higher density 
molecular g 10 4 /r cm 3 ) and adopted a HCN-to-H 2 conversion factor 

^dense(-f^2)/-^HCN ~ 10 M Q (K kms" 1 pc 2 )^ 1 . Luminous and ultraluminous infrared 
galaxies usually have much higher HCN brightness temperature, which might result in 
a lower conversion factor. Analyses including LVG calculations indicate that the HCN 
luminosity measures the mass of dense molecular gas and the HCN/CO luminosity ratio 
indicates the fraction of molecular gas in a dense phase. 

We thank Simon Radford, Dennis Downes, and Mark Heyer for help with some of the 
observations. We appreciate the generous support and allocation of observing time from 
the NRAO 12m, the IRAM 30m, and the FCRAO 14m. We would also like to acknowledge 
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review and helpful comments of the anonymous referee are greatly acknowledged. This 
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35.5 


2.8 


20 


27 


SAB(rs)c 


Arp 55 


09 12 38.1 


44 32 26 


11957 


162.7 


0.7 






GPair,HII 


NGC 2903 


09 29 19.9 


21 43 19 


556 


6.2 


12.6 


17 


65 


SB(s)d,HII 


*UGC 05101 


09 32 04.8 


61 34 37 


11810 


160.2 


1.1 






S?LINER,Syl.5 


M82 


09 51 42.9 


69 54 59 


203 


3.5 


11.2 


65 


66 


I0,Sbrst,HII 


NGC 3079 


09 58 35.0 


55 55 17 


1150 


16.2 


7.9 


165 


80 


SB(s)c,LINER,Sy2 


*IRAS 10565+2448 


10 56 36.2 


24 48 40 


12923 


173.3 


0.4 






LINER, HII 



Table 1. (continued) 



Galaxies a 


R.A. (1950) 


DEC. 


yfkms- 1 )' 1 


d(Mpc) c 


D-25 (') 


P.A. d 


Inclination 


Classifications 6 


Arp ±4o 


11 f|1 n p: e 

11 Ul uo.o 


a 1 07 1 n 

41 U/ 1U 


iuoy4 


1 A 1 1 
140.0 


U.O 






Ring,GPair 


Ly\j\^ ODDO 


11 Uo oO.o 


DO DO OO 


oyo 


in c 
1U.0 


Q Q 
o.o 


ou 


7fi 


OJJ^S )CQ 


WPP 1£97 
l^'j^ 00Z 1 


11 17 17 O 

ii i / tS / .y 


11 1C HQ 
±0 ID Uo 


oyo 


7 £. 


O 1 

y.i 


1 71 
I/O 


DO 


DAISYS ) D,oyz 


NGC 3628 


11 1 7 fi 


13 51 48 


847 


9.6 


14.8 


104 


78 


sti "npc"sn 


NGC 3893 


11 46 01.1 


48 59 20 


968 


13.9 


4.5 


165 


55 


SAB(rs)c 


NGC 4030 


11 57 50.3 


-00 49 22 


1460 


17.1 


4.2 


30 


45 


SA(s)bc 


NGC 4041 


11 59 39.0 


62 24 56 


1234 


18.0 


2.7 




21 


SA(rs)bc 


NGC 4414 


12 23 57.9 


31 30 00 


720 


9.3 


3.6 


155 


56 


SA(rs)c? 


NGC 4631 


12 39 41.5 


32 48 54 


620 


8.1 


15.5 


86 


80 


SB(s)d 


NGC 4826 


12 54 16.9 


21 57 18 


408 


4.7 


10.0 


115 


57 


(R)SA(rs)ab,Sy2 


NGC 5005 


13 08 38.7 


37 19 23 


950 


14.0 


5.8 


65 


61 


SAB(rs)bc,Sy2 


NGC 5055 


13 13 36.4 


42 17 34 


504 


7.3 


12.6 


105 


55 


SA(rs)bc,HII,LINER 


NGC 5135 


13 22 56.5 


-29 34 26 


4114 


51.7 


2.6 




45 


SB(l)ab,Sy2 


M83 


13 34 11.0 


-29 36 48 


518 


3.7 


12.9 


45 


27 


SAB(s)c,HII,Sbrst 


*Mrk 273 


13 42 51.7 


56 08 14 


11324 


152.4 


1.1 






Sy2, LINER 


NGC 5678 


14 30 37.1 


58 08 35 


1922 


27.8 


3.3 


5 


61 


SAB(rs)b 


NGC 5713 


14 37 37.5 


-00 04 35 


2180 


24.0 


2.8 


10 


27 


SAB(rs)bc:pec 


NGC 5775 


14 51 26.8 


03 44 53 


1690 


21.3 


4.2 


146 


78 


SBc?:sp 


*IRAS 17208-0014 


17 20 47.8 


-00 14 15 


12835 


168.0 








HII,Sbrst 


IRAS 18293-3413 


18 29 21.4 


-34 13 42 


5449 


72.1 








HII 


NGC 6701 


18 42 35.5 


60 36 08 


3950 


56.8 


1.5 


19 




(R')SB(s)a,HII 


NGC 6921 f 


20 26 20.8 


25 33 23 


4399 




0.9 


141 


74 


SA(r)0/a 



Table 1. (continued) 



Galaxies 8 


R.A. (1950) 


DEC. 


V(kms- 1 ) h 


d(Mpc) c 


£>25 (') 


P.A. d 


Inclination 


Classifications' 5 


NGC 6946 


20 33 49.5 


59 58 47 


48 


5.5 


11.5 


19 


30 


SAB(rs)cd,HII 


NGC 7130 


21 45 19.9 


-35 11 01 


4842 


65.0 


1.5 






Sa:pec,Sy2,LINER 


IC 5179 


22 13 13.0 


-37 05 30 


3447 


46.2 


2.3 


57 


61 


SA(rs)bc,HII 


NGC 7331 


22 34 46.7 


34 09 21 


821 


15.0 


10.5 


90 


69 


SA(s)b,LINER 


NGC 7469 


23 00 44.4 


08 36 16 


4812 


67.5 


1.5 


125 


46 


(R')SAB(rs)a,Syl.2 


NGC 7479 


23 02 26.8 


12 03 06 


2385 


35.2 


4.1 


25 


39 


SB(s)c,LINER,Sy2 


*IRAS 23365+3604 


23 36 32.3 


36 04 33 


19330 


266.1 


0.5 






Liner 


Mrk 331 


23 48 53.5 


20 18 27 


5363 


75.3 


0.7 






S?,HII,Sy2 



Notes to Table 1. 

Galaxy names marked with * indicates ULIGs with Lir^10 12 L q . Galaxies with HCN data available in the literature that have been included in our 
HCN sample are not listed (see Table 1 in the companion paper [Gao & Solomon 2003, Paper II] for these galaxies). 

a Galaxies that were mapped along the major axes are marked in boldface and the details will be discussed in a future paper (Gao & Solomon in 
preparation). 

b Heliocentric velocity drawn from the literature and NED or determined from our CO observations. 

c Source distance is calculated using Ho — 75 kms" 1 Mpc -1 corrected for the local group motion, except for a few nearby galaxies, where the most 
recent and commonly used values in the literature are adopted. 

d Position angles of the major axes, except for NGC 1068 and IC 342, that are used in mapping HCN emission. 
e Almost entirely from NED, GPair— galaxy pair, GGroup— galaxy group, Sbrst— starburst galaxy, Ring— ring galaxy. 

f NGC 6921 is most likely only partially associated with the IRAS source IRAS 20264+2533, which is possibly more related to the galaxy MCG 
+04-48-002, - 1'.5 east of NGC 6921 (Sanders et al. 2003). NGC 6921 was included in our HCN survey since Young et al. (1986) reported relatively 
strong CO emission and assumed that the IRAS emission comes entirely from NGC 6921. 
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Table 2. Telescope Efficiencies at 3 mm a 

Telescope Typical T sys T Scale ^ b F cS c B cS c r]B d 

NRAO 12m 220 (250)K T£ 0.89 (0.83) 

IRAM 30m 250 (270)K T* A 0.92 (0.92) 0.75 (0.65) 

FCRAO 14m 500 (850)K T* A 0.60 (0.55) 



a Both values at 3.4 mm for HCN J=l-0 and at 2.6 mm for CO J=l-0 (in parentheses) are 
given. These efficiencies were adopted for conversion from the measured antenna temperature 
scale to the main-beam T m b scale. 

b 2mb = Tr/^m; where rj^ is the corrected main-beam efficiency. 

c Tmb = T A (F e s/B eS ), where F eS and B cf[ are the forward and effective beam efficiencies 
respectively. 

d T m i, = T^/t/b, where r/ B = B cS /F cS is the main-beam efficiency. 



Table 4. Observed Quantities in the HCN Survey 



Galaxies 


Offsets 
a 


Ihcn 


"HCN 

K km 


Ico 
s- 1 


"CO 


Vco 
km s" 


Ay 

-i 




/lOO^m 

Jy 


J ^HCN^i 

£3 mb 

Kkm 


£2 mb 

s- 1 


ScoAF ShcnAI/ 
Jy km s _1 


NGC 253 


map a,b 


23.0 


0.6 


272.0 


1.6 


251.9 


202.4 


0.54 


1857.8 


33.5 


850.5 


24523.0 


947.5 




map a,c 


74.2 


1.5 














153.8 






770.5 


IC 1623 b 


0, 


.81 


.15 


20.6 


0.4 


6076 


103 


0.77 


30.6 


.91 


24.8 


745.0 


27.3 


NGC 660 b 


0, 


.85 


.20 


31.7 


1.6 


825 


316 


0.63 


103.8 


>0.96 


45.0 


>1145.8 


>28.0 


NGC 695 c 


0, 


1.4 


.24 






9727 


165 


0.55 


14.2 


1.4 


7.3 d 


214.5 


6.7 


Mrk 1027 c 


0, 


0.91 


.17 


29.1 


.90 


9007 


150 


0.64 


8.5 


.91 


29.1 


136.0 


4.4 


NGC 891 


map a,e 


0.85 


0.42 


35.5 


2.1 


548 


160 


0.31 


198.6 


<2.71 


158.8 


3733.7 


<63.6 




map a,c 


4.8 


0.5 














14.5 






71.5 


NGC 1022 b 


0, 


.28 


.09 


9.6 


.37 


1463 


89 


0.75 


26.8 


0.32 


11.6 


337.7 


9.4 


NGC 1055 b 


0, 


<.96 


.48 










0.36 


63.8 


<1.1 


75.3 d 


2221.4 


<31.3 


NGC 1068 


map a,e 


6.6 


0.6 


69.0 


2.8 


1185 


175 


0.79 


251.8 


13.0 


155.7 


3651.1 


302.2 


NGC 1144 c 


0, 


.82 


0.24 


46.2 


1.0 


8750 


650 


0.46 


11.6 


0.82 


46.2 


216.9 


3.5 


NGC 1365 


map a,b 


2.2 


0.3 


55.7 


1.6 


1623.4 


233.4 


0.45 


185.4 


5.5 


168.0 


3710.0 


179.8 


IC 342 


map a,b 


3.4 


0.3 


50.0 


0.8 


35.2 


72.8 


0.39 


660.7 


18.0 


808.2 


18962.0 


522.0 




map a,e 


3.2 


0.3 


49.7 


0.8 


32.0 


70.0 














NGC 1614 c 


0, 


1.5 


.22 


43.5 


1.2 


4750 


184 


0.96 


33.0 


1.5 


43.5 


208.8 


7.2 


VII Zw 31 c 


0, 


1.1 


.17 


21.0 


0.9 


16275 


200 


0.55 


9.8 


1.1 


21.0 


99.0 


5.3 


IRAS 05189-2524° 


0, 


1.1 


.20 






12831 


110 


1.1 


12.0 


1.1 


3.1 d 


92.4 


5.3 


NGC 2146 


map a,e 


1.9 


0.3 


50.1 


1.8 


891 


311 


0.71 


204.2 


<4.5 


105.7 


2491.2 


<106.2 




map a,c 


8.1 


0.6 














21.2 


413.2 


2400.0 


104.1 


NGC 2276 b 


0, 


.23 


.06 


8.3 


.19 


2405 


61 


0.42 


30.8 


0.25 


10.1 


292.6 


6.8 


Arp 55 c 


0, 


.81 


.15 






11945 


165 


0.59 


10.1 


0.81 


6.1 d 


181.3 


3.8 


NGC 2903 b 


0, 


1.7 


0.25 


27.8 


0.8 


564.5 


155.6 


0.36 


147.2 


>1.7 


~ 72.0 


2105.0 


48.6 




map a,c 


4.4 


0.5 














~ 10.0 


~ 414.0 d 


1987.2 


48.0 


UGC 5101 c 


0, 


2.1 


.33 


15.6 


0.7 


11785 


245 


0.59 


20.4 


2.1 


15.6 


73.3 


10.4 


M82 


map a,b 


6.6 


0.4 


192.7 


2.1 


235.5 


180.5 


0.94 


1349.0 


15.6 


607.7 


17671.9 


458.7 




map a,c 


27.4 


1.4 














93.0 






445.2 




map a ' e 


8.2 


0.7 


257.6 


1.9 


198.5 








22.0 


660.0 


15648.5 


507.2 



Table 4. (continued) 



( -i 1 Yl £»Q 
V_T Cll tX A. 1 C o 


Wllbt Lo 
II 


JHCN 


^HCN 

K km s 


^co 

-l 


CT CO 


km s 


Ay 

-i 


/lOO^m 


JWOftm 

Jy 


Kkm s 


f IcodQ, 
-1 


ScoAV ShcnAV 
Jy km s _1 


NGC 3079 b 


0, 


2.6 


.42 


66.1 


1.1 


1111 


365 


0.48 


103.5 


2.9 


110.0 


3187 5 


77 6 


IRAS 10565+2448° 


0, 


1.8 


.20 


16.5 


0.6 


12950 


380 


0.85 


14.4 


1.8 


16.5 


77 ^ 
i i .0 


O.D 


Arp 148 c 


0, 


1.1 


.16 






10520 


220 


0.57 


10.7 


1.1 


18.3 d 


86.0 


5.4 


NGC 3556 b 


0, 


.50 


.08 


10.1 


.80 


674 


89 


0.40 


80.7 


.57 


40.0 


1032.5 


16.8 


NGC 3627 b 


0, 


.68 


.12 






687 


159 


0.39 


144.9 


>0.76 


190. 7 d 


4477.4 


>21.4 


NGC 3628 


map a,b 


1.2 


0.2 


57.8 


1.1 


844.1 


193.7 


0.40 


121.9 


3.1 


185.5 


4351.8 


90.2 


NGC 3893 b 


0, 


.80 


.13 


8.9 


.34 


968 


152 


.39 


38.2 


0.90 


26.6 


626.5 


26.0 


NGC 4030 b 


0, 


1.2 


.26 


20.4 


.68 


1458 


232 


0.38 


50.4 


1.35 


65.1 


1879.9 


39.3 


NGC 4041 b 


0, 


.35 


.07 


12.2 


.48 


1212 


116 


0.43 


32.6 


0.40 


15.0 


433.6 


11.7 


NGC 3079 b 


0, 


2.6 


.42 


66.1 


1.1 


1111 


365 


0.48 


103.5 


2.9 


110.0 


3187.5 


77.6 


NGC 4414 b 


0, 


1.2 


0.3 


30.4 


.70 


706 


269 


0.38 


68.4 


-1.0 


66.3 


1982.0 


38.5 


NGC 4631 b 


0, 


.72 


.16 


12.1 


.95 


658 


75 


0.40 


208.4 


>0.92 


40.5 


1167.0 


>27.1 




map a ' e 


1.1 


0.3 


18.5 


1.2 


603 


125 






2.1 


49.6 


1166.0 


49.3 


NGC 4826 b 


0, 


1.0 


.17 


22.2 


1.6 


372 


216 


0.49 


75.7 


>1.4 


75.0 


2172.6 


>40.1 


NGC 5005 b 


0, 


1.3 


.25 


18.7 


1.5 


1034 


210 


0.32 


59.4 


2.6 


87.7 


2542.3 


75.4 




45, 21 


0.75 


.25 


10.0 


1.5 


1176 


60 














NGC 5055 b 


0, 


1.5 


.15 


18.8 


1.3 


500 


226 


0.25 


157.4 


>1.6 


200.1 


4690.5 


>46.5 


NGC 5135 b 


0, 


.70 


.08 


11.0 


.45 


4111 


92 


0.58 


29.0 


0.74 


14.5 


426.3 


21.2 


M83 


map a,b 


63.7 


5.4 


2.5 


0.3 


508.5 


75.1 


0.42 


638.3 


12.9 


557.1 


13085.2 


376.6 


Mrk 273 c 


0, 


3.5 


.47 


22.7 


.89 


11320 


357 


1.1 


22.3 


3.5 


22.7 


107.0 


16.8 


NGC 5678 b 


0, 


.71 


.22 


10.8 


1.0 


1928 


266 


0.36 


26.7 


.81 


28.0 


811.7 


23.1 


NGC 5713 b 


0, 


.25 


.07 


14.4 


.30 


1913 


97 


0.28 


37.2 


0.48 


17.5 


505.9 


10.0 


NGC 5775 b 


0, 


.80 


.25 


12.5 


1.4 


1697 


151 


0.46 


51.4 


~0.90 


35.6 


868.7 


26.5 



Table 4. (continued) 



Galaxies 


Offsets 

II 


^HCN 


"HCN 

K kms 


-fco 

-l 


"CO 


Vco 
km s 


Ay 

-i 


/lOOpm 


/lOOfim 

Jy 


(' T WO 

S2 mb 

Kkms 


J IcodCls 

Simb 

-1 


ScoAV ShcnAV 
Jy km s _1 


IRAS 17208-0014 b 


0, 


.91 


.19 


5.7 


0.3 


12855 


303 


0.98 


34.9 


1.0 


6.9 


206.5 


30.5 


IRAS 18293-3413 b 


0, 


.50 


.16 


16.8 


1.0 


5490 


210 


0.69 


53.8 


0.56 


20.4 


591.8 


16.1 


NGC 6701 b 


0, 


.28 


.06 


10.8 


.20 


3946 


88 


0.48 


21.4 


0.31 


13.1 


380.9 


9.0 


NGC 6921 b 


0, 


.50 


0.1 


5.0 d 


.31 


4200 


397 


0.54 


18.5 


0.56 


6.0 


182.0 


16.5 


NGC 6946 


map a,b 


2.6 


0.2 


58.5 


1.4 


37.7 


125.0 


0.40 


344.4 


11.5 


380.0 


11400.5 


327.8 




map a,c 


11.4 


1.0 














59.5 






292.6 


NGC 7130 b 


0, 


.50 


.11 


10.9 


.90 


4842 


87 


0.64 


26.6 


0.57 


13.2 


383.0 


16.5 


IC 5179 b 


0, 


1.1 


.20 


12.7 


1.6 


3418 


242 


0.50 


38.5 


1.2 


15.4 


446.6 


37.1 


NGC 7331 b 


0, 


1.3 


.27 


16.0 


1.1 


871 


198 


0.31 


114.8 


>1.5 


>59.4 


>1752.3 


<42.8 


NGC 7469 b 


0, 


.31 


.06 


8.3 


.CO 


4925 


177 


0.78 


35.4 


0.35 


10.1 


293.8 


10.5 


NGC 7479 b 


0, 


.60 


.12 


11.0 


.54 


2368 


192 


0.59 


26.2 


0.67 


26.7 


787.6 


20.2 


IRAS 23365+3604° 


0, 


1.2 


.24 


10.9 


0.5 


19300 


68 


0.94 


8.2 


1.2 


10.9 


51.2 


5.9 


Mrk 331 b 


0, 


.38 


.11 


8.6 


.65 


5410 


215 


0.81 


21.0 


0.43 


11.4 


336.3 


12.9 



Notes to Table 4. 

The parameters/quantities listed are: galaxy name (column 1); observed position (off-sets or major-axis mapping) as given in Table 1 (column 2); measured 
line intensities and estimated errors for HCN and CO respectively in columns 3 to 6, but CO data taken from literature are not listed; heliocentric velocity 
of the galaxy as determined from either our CO observations or CO data in literature and the linewidth (FWHM) in column 7, 8; the IRAS far-IR flux 
density ratio jeo^m/Aoo^m an d 100 fim flux (column 9, 10); the integrated line intensities over the entire source on the main-beam antenna temperature 
scale for HCN and CO (column 11, 12); and the total CO and HCN line emission fluxes (last two columns). 

a Multiple position HCN and CO observations along the major axes are available (Gao & Solomon in preparation). Here only the central position (0, 
0) data are given. 

b Data from the NRAO 12m are on the temperature scale of T^. Only NGC 1055 was not detected and is listed as a 2a upper limit. However, there 
are still limits in the total integrated HCN fluxes in some nearby galaxies because of either only one measurement towards the center (lower limits) or 
non-detections at large galactic radii but detected at the centers/inner disks (upper limits). 

c Data from the IRAM 30m are converted to the main beam temperature scale. 

d CO data from literature: NGC 695, IRAS05189-2524, Arp 55 from Sanders et al. (1991); Arp 148 from Horellou et al. (1995); NGC 1055 from Sage 
(1987); NGC 1365 from Sandqvist, Jorsater, & Lindblad (1995); NGC 2903 from Jackson et al. (1991); NGC 6921 from this work, but Young et al. (1986, 
1995) listed the line intensity twice as large as ours (see more notes for NGC 6921 in Table 1). 

e Data from the FCRAO 14m are expressed on the main beam temperature scale. In most cases, only the central beams covering the nuclear regions 
have HCN detections. The total integrated HCN line fluxes have been derived using the weak detections or the 2a upper limits from other beams covered 
by QUARRY along the major axes along with the central beam measurements. 



Table 4. Global Properties of Galaxies in the HCN Survey 



Galaxies 




Leo 


T a 
-^HCN 


LhCN/ LnO 


Lib b 


Lib 


M dcnBC (H 2 ) 


M(H 2 ) 




10 10 L n 

-LW -LJQ 


10 8 K kins" 1 pc 2 




L n /Mn 
J - / ©/ lvJ -0 




NGC 253 


2.1 


4.6 


0.27 


0.059 


77.8 


9.6 


IC 1623 


46.7 


130.5 


8.5 


0.065 


55.0 


7.5 


NGC 660 


3.7 


7.3 


>0.26 


>0.036 


<142.3 


10.6 


NGC 695 


46.6 


92.9 


4.3 


0.046 


108.4 


10.5 


MRK 1027 


25.7 


41.7 


1.89 


0.045 


136.0 


12.9 


NGC 891 


2.6 


11.0 


0.25 


0.024 


104.0 


5.4 


NGC 1022 


2.6 


4.2 


0.20 


0.047 


131.5 


13.1 


NGC 1055 


2.1 


13.3 


<0.37 


<0.028 


>56.9 


3.3 


NGC 1068 


28.3 


20.7 


3.61 


0.174 


78.5 


28.6 


NGC 1144 


25.1 


108.9 


2.67 


0.025 


94.0 


4.8 


NGC 1365 


12.9 


58.7 


3.10 


0.053 


41.6 


4.6 


IC 342 


1.4 


9.4 


0.47 


0.050 


30.2 


3.4 


NGC 1614 


38.6 


24.5 


1.25 


0.051 


306.8 


33.2 


*VIIZw31 


87.1 


125.0 


9.8 


0.078 


88.9 


14.6 


*05189-2524 


118.1 


67.0 


6.2 


0.093 


190.4 


36.9 


NGC 2146 


10.0 


12.5 


0.96 


0.071 


122.4 


16.8 


NGC 2276 


6.2 


10.2 


0.40 


0.039 


156.0 


12.8 


ARP 55 


45.7 


125.0 


3.8 


0.030 


120.2 


7.6 


NGC 2903 


0.83 


2.3 


>0.09 


>0.036 


<104.0 


7.6 


*UGC 05101 


89.2 


50.8 


10.0 


0.197 


89.2 


36.7 


M82 


4.6 


5.7 


0.30 


0.053 


153.3 


16.9 


NGC 3079 


4.3 


24.0 


~1.0 


~0.042 


-43.5 


3.7 


*10566+2448 


93.8 


61.5 


10.2 


0.166 


92.0 


32.0 


ARP 148 


36.5 


>47.0 


4.0 


<0.085 


91.2 


<16.2 


NGC 3556 


1.35 


>4.5 


>0.09 


-0.020 


<155.0 


<6.3 


NGC 3627 


1.26 


4.4 


>0.08 


>0.017 


<160.0 


6.0 


NGC 3628 


1.01 


7.1 


0.24 


0.034 


43.5 


3.3 


NGC 3893 


1.15 


4.1 


0.23 


0.056 


50.4 


5.9 


NGC 4030 


2.14 


15.2 


0.54 


0.036 


40.0 


3.0 


NGC 4041 


1.70 


3.9 


0.18 


0.046 


94.8 


9.1 


NGC 4414 


0.81 


4.6 


0.16 


0.033 


54.8 


3.7 


NGC 4631 


2.0 


2.3 


-0.09 


~0.039 


-222.2 


18.5 


NGC 4826 


0.26 


1.3 


>0.04 


>0.030 


<65.5 


4.1 


NGC 5005 


1.4 


8.2 


0.41 


0.049 


21.5 


2.2 
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Table 4 — Continued 



Galaxies 






T a 
-^HCN 


Lrcn/ Leo 


Lin b 

""dense V ^ I 


M(H 2 ) 




10 10 L 


10 8 K kms- 1 pc 2 




L /M 




NGC 5055 


1.1 


8.6 


>0.10 


>0.012 


<111.0 


2.7 


NGC 5135 


13.8 


31.3 


2.73 


0.087 


50.6 


9.1 


M83 


1.4 


8.1 


0.35 


0.043 


40.5 


3.7 


*MRK 273 


129.9 


65.0 


15.2 


0.234 


85.5 


41.8 


NGC 5678 


3.0 


17.2 


0.75 


0.044 


40.0 


3.6 


NGC 5713 


4.2 


8.1 


0.22 


0.027 


188.9 


10.7 


NGC 5775 


3.8 


10.9 


0.57 


0.052 


66.7 


7.2 


*17208-0014 


234.5 


146.9 


37.6 


0.256 


62.4 


33.4 


18293-3413 


53.7 


85.5 


4.03 


0.047 


66.6 


13.2 


NGC 6701 


11.2 


34.0 


1.38 


0.041 


40.6 


6.9 


NGC 6921 


11.4 


17.5 


~2.81 


~0.160 


~41.0 


13.6 


NGC 6946 


1.6 


9.2 


0.49 


0.053 


33.5 


3.6 


NGC 7130 


21.4 


44.9 


3.27 


0.071 


65.4 


10.0 


IC 5179 


14.1 


-26.4 


3.42 


~0.129 


41.2 


~11.2 


NGC 7331 


3.5 


>10.7 


>0.44 


~0.041 


<79.8 


<6.8 


NGC 7469 


40.7 


37.1 


2.19 


0.059 


185.4 


23.0 


NGC 7479 


7.4 


26.7 


1.12 


0.042 


66.1 


5.8 


*23365+3604 


142.0 


85.0 


15.0 


0.176 


94.6 


35.0 


MRK 331 


26.9 


52.1 


3.35 


0.064 


80.3 


10.7 



a The 2 a upper limit is listed for NGC 1055. The lower limits are for nearby galaxies 
where we only detected HCN in the galaxy centers or extensive mapping is still required in 
order to estimate the total HCN and CO luminosities. 

b M densc = IOLhcn M /K kms" 1 pc 2 and M(H 2 ) = 4.78L co M /K kms" 1 pc 2 , see §4.5. 
Note. — ULIGs with L IR ^,1O 12 L are indicated by asterisks (*). 
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Fig. 1. — (a) The distribution of infrared luminosities in our HCN survey sample of 53 
galaxies plus a dozen galaxies with published HCN observations. More than half of the 
sample galaxies have Ljr < 10 11 L Q and most of them have < 3 x 1O 1O L . (b) The 
distribution of CO luminosities in the same 65 galaxies, (c) The distribution of HCN 
luminosities in the same sample. 

Fig. 2. — The HCN spectra observed with the IRAM 30m telescope. The vertical axis is the 
main-beam temperature scale in units of mK. For nearby galaxies, only the central spectra 
are shown here. CO spectra (dashed lines) from the same telescope have been divided by 
indicated factors for comparison. The spectra are ordered by right ascension, as in Table 1. 
Some spectra are CO J=2-l, which have a much smaller beam size than the HCN J=l-0 
observations. Some frames have only the CO line emission windows (derived from published 
data) indicated. For NGC 253, we show the CO J=2-l spectrum obtained from the NRAO 
12m (scaled to the main-beam temperature scale for proper comparison), which has roughly 
the same beam size as that of the HCN J=l-0 from the IRAM 30m. 

Fig. 3. — HCN spectra obtained from the NRAO 12m telescope. The temperature scale is 
Tjj in mK. T m & = T^/rj* M , with a corrected beam efficiency rf M ~ 0.89. As in Fig. 2, only the 
central spectra are shown here for nearby galaxies and a few of them (NGC 253, NGC 2903, 
M82, and NGC 6946) are the observations done at both the 12m and 30m telescopes. The 
dashed lines are CO J=l-0 spectra divided by indicated factors. The ultraluminous infrared 
galaxy IRAS 17208—0014 (the only one observed with the 12m) is shown last. The spectra 
of NGC 1068 (obtained from the FCRAO 14m) and NGC 1055 (not detected) are not shown. 

Fig. 4. — HCN and IR luminosities in our HCN survey sample {open circles) compared 
with the original correlation proposed by Solomon et al. (1992) for only 10 galaxies (circles 
with asterisks). Note that we used the IR luminosity rather than the far-IR luminosity as in 
Solomon et al. (1992). Representative error bar is shown for the highest Lhcn in the sample, 
which includes ~20% calibration uncertainty. The line indicates logLi R = 0.97/ogL H CN + 3.1 
(see text). 

Fig. 5. — (a) Relationship between the HCN luminosity and the luminosity ratio of the 
12 /im-to-IR emission, (b) The HCN-to-CO luminosity ratio is plotted against the ratio of 
the 12-to-100/xm emission. There are no significant correlations found implying that galaxies 
with excess mid-IR emission do not show stronger HCN emission and do not contain higher 
fractions of dense molecular gas. The sample is divided into luminous and ultraluminous 
galaxies (LIGs/ULIGs) with LpR > 10 11 L o (filled circles) and less luminous galaxies (open 
circles). Representative error bar is also shown for the highest Lhcn- 
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Fig. 2a See 0310341.fg2a.gif 



Fig. 2b 

See 0310341.fg2b.gif 



Fig. 2c 

See 0310341.fg2c.gif 



Fig. 3a 

See 0310341.fg3a.gif 



Fig. 3b 

See 0310341.fg3b.gif 



- 36 - 




-37- 




- 38 - 




- 39 - 




. 1 1 1 1 1 1 1 1 1 1 1 

" 1720-00 f Dn 

r n r 

: Pi 'i 

r Ml 

J 1 L 

: I" 1 J 
- 1 " i , 


i i 

CO/5 

H 


i : 


l~Lt L 
■ i i i i i i i i i i i 




i 
i : 



1.25 10 4 1.3 10 4 1.35 10 4 

Velocity (km/s) 



-41 - 





-2.5 -2.0 -1.5 -1.0 -0.5 0.0 

lo g( L 12um / 



0.0 



~i 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



~i 1 1 1 1 r 



LIGs/ULIGs L IR >10 U L 







O Normal Spirals L IR <10 11 L Q 



0.5 



1.0 



1.5 



2.0 



2.5 



• • • 

o oo ojgjt 




J I I I I I I I I I I I I I L_ 



_l I I I I L_ 



-3.0 -2.5 -2.0 -1.5 -1.0 -0. 

lo g( f 12 M m/ f 100 /zm / 



This figure "fg2a.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0310341v2 



This figure "fg2b.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0310341v2 



This figure "fg2c.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0310341v2 



This figure "fg3a.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0310341v2 



This figure "fg3b.gif" is available in "gif" format from: 



http://arXiv.org/ps/astro-ph/0310341v2 



